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> ■ Abstract 
in 

Differential cross sections and complete sets of polarization observables are presented for the 
Gamow- Teller 12 C(p, n) 12 N(g.s., 1 + ) reaction at a bombarding energy of 296 MeV with momentum 
transfers q of 0.1 to 2.2 fm . The polarization transfer observables are used to deduce the spin- 
longitudinal cross section, ID q , and spin-transverse cross sections, ID P and ID n . The data are 



compared with calculations based on the distorted wave impulse approximation (DWIA) using 



shell-model wave functions. Significant differences between the experimental and theoretical results 
are observed for all three spin-dependent IDi at momentum transfers of q > 0.5 fm -1 , suggesting 
the existence of nuclear correlations beyond the shell model. We also performed DWIA calculations 
employing random phase approximation (RPA) response functions and found that the observed 
discrepancy is partly resolved by the pionic and rho-mesonic correlation effects. 
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I. INTRODUCTION 



Nuclear spin-isospin correlations produce various interesting phenomena in nuclei de- 
pending on the momentum transfer q jl|. At small momentum transfers, q ~ fm _1 , the 
spin-isospin particle-hole interaction is strongly repulsive, which leads to such collective ex- 
citations in nuclei as the Gamow- Teller (GT) giant resonance. The quenching of the total 
strength of the GT transitions from its sum rule value has prompted theoretical studies of 
possible mechanisms, ranging from conventional configuration mixing to admixture of the 



A-hole (A-h) states. Recent experimental studies of Zr(p,n) and (n,p) reactions 
have revealed that the coupling to 2-particle-2-hole (2p-2h) excitations is the major source 
of quenching of the GT strengths, while A-h coupling plays a minor role. 

At fairly large momentum transfers, q > 1 fm _1 , the spin-transverse interaction induced 
by one rho-meson exchange remains repulsive, while the spin-longitudinal interaction driven 
by one pion exchange becomes attractive. This attraction of the spin-longitudinal inter- 
action produces pion condensation and its precursor phenomena. Pion condensation 4jlis 

n i 

expected to occur in neutron stars (NSs) such as 3C58 |5j and accelerate their cooling |6j. 
It is predicted that pion condensation does not occur in normal nuclei. However, precursor 
phenomena may be observed in normal nuclei if they are near the critical point of the phase 
transition. As possible evidence of a precursor, enhancements of the Ml cross section in pro- 



ton inelastic scattering [3,0] and of the ratio Rl/ Rt, the spin-longitudinal (pionic) response 
function Rl to the spin-transverse (rho-mesonic) response function Rt, in the quasielastic 
scattering (QES) region were proposed at a momentum transfer of about q ~ 1.7 fm . 

Motivated by these predictions, many experiments involving the Ml transition and QES 



have been carried out. These include the measurement o 



the 12 Cip,p') 12 C*(l+T = i) 



scattering at incident energies of about 120 to 800 MeV ll|, [12|, [13J, ll4| , and the experimen- 
tal extraction of RtJRt u sing complete sets of polarization transfer observables in (p,p ') 



scattering [15 



16 



17 



18 



19] and (p, n) reactions [20j, 



21 



22 



23 



24j on various targets 



at intermediate energies. However, these experimental data did not show the theoretically 
expected enhancements. Recent analysis of QES data shows pionic enhancement in the spin- 
longitudinal cross section representing Rl, which suggests that the lack of enhancement in 
the value of Rl/Rt is due to the rho-mesonic component [25j. The measurement of the pure 
pionic excitation of 16 O(p,p') 16 O*(0~, T = 1) scattering at T p = 295 MeV also supports such 
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an enhancement 



26|. 



The analyses of QES and pure pionic excitation data suggest that one needs to recon- 
sider the interpretation of the Ml data. The pionic effect in the Ml cross section might 
be masked by the contribution from the rho-mesonic component because the Ml state is a 
mixture of spin-longitudinal and spin-transverse states. In addition, proton inelastic scat- 
tering, including the Ml transition, mix isoscalar and isovector contributions. With respect 
to this issue, the GT 12 C(p, n) 12 N(g.s., 1 + ) reaction, which is the isobaric analog of the Ml 
12 C(p,p') 12 C*(l + , T = 1) scattering, has an advantage because the (p, n) reaction exclusively 
depends on the isovector contribution. Thus, it is interesting to study nuclear correlation 
effects in this reaction by separating the cross section into pionic and rho-mesonic compo- 
nents, using a complete set of polarization observables. However, no complete polarization 
transfer measurements for the GT reaction have been reported until now. 

In this paper, we present differential cross sections and complete sets of polarization 
observables for the 12 C(p, n) 12 N(g.s., 1 + ) reaction measured at T p = 296 MeV over a range 
of q = 0.1 to 2.2 fm -1 . This incident energy is ideal for studying GT transitions because 
the spin excitations are dominant in the (p,n) reaction near 300 MeV {27]. In addition, 
distortion effects are minimal around 300 MeV [1], thereby enabling one to extract reliable 
nuclear structure information such as nuclear correlation effects. 

The differential cross section and a complete set of polarization transfer observables are 
used to separate the cross section into nonspin (I Do), spin-longitudinal (ID q ), and two 
spin-transverse (ID P and ID n ) polarized cross sections. The spin-dependent polarized cross 
sections, ID q , ID P , and ID n , are compared with distorted wave impulse approximation 
(DWIA) calculations using random phase approximation (RPA) response functions, in order 
to assess the nuclear correlation effects quantitatively. The theoretical calculations give good 
descriptions of all of the spin-dependent polarized cross sections. These results demonstrate 
the existence of pionic and rho-mesonic correlations inside nuclei. 



II. EXPERIMENTAL METHODS 



The experiment was performed using the Neutron Time-Of-Flight (NTOF) facility 28] 
and the neutron detector/polarimeter NPOL3 [29] at the Research Center for Nuclear 
Physics (RCNP) at Osaka University. Detailed descriptions of the NTOF facility and the 
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NP0L3 system can be found in Refs. 



28 



29 



30 



31 



32 



33 



34 



35] . Only the details relevant 



to the present experiment are described here. Schematic layouts of the NTOF facility and 
the NPOL3 system are presented in Figs. [T] and [2J respectively. In Fig. [TJ the coordinate 
systems for incident protons and outgoing neutrons are also shown. 



A. Polarized proton beam 

The high intensity polarized ion source (HIPIS) [36] at RCNP was used to provide the 
proton beam. The direction of the beam polarization was reversed every 5 s in order to 
minimize geometrical false asymmetries that might be present in the experimental apparatus. 
The polarized proton beam from HIPIS was accelerated up to T p = 53 and 296 MeV by 



38 



cyclotrons, respectively. The radio 



the Azimuthally Varying Field (AVF) [37| and Ring 
frequency (rf) of the AVF cyclotron was 15.42 MHz, corresponding to a beam pulse period 
of 64.86 ns. In the cross section and analyzing power measurement, one out of every five 
beam pulses was selected prior to injection into the Ring cyclotron, yielding a pulse period 
of 324.3 ns. This pulse selection reduced the wrap-around events of slow neutrons from 
preceding beam pulses. In the measurement of the polarization transfer observables, the 
pulse selection was not performed, so as to achieve reasonable statistical accuracy. Note 
that the contribution from the wrap-around events was negligibly small (< 1%) because of 
the kinematical selection at the NPOL3 system (see Sec. IIII Cl) . The single-turn extraction 
was maintained during these measurements in order to preserve the beam polarization. 

Two superconducting solenoid magnets (SOLI and SOL2) located in the injection line 
from the AVF to Ring cyclotrons were used to precess the proton spin. Each magnet 
can rotate the direction of the polarization vector from the normal N into the sideways S 
directions. These two magnets are separated by a bending angle of 45°, and thus the spin 
precession angle in this bending magnet is about 85.2° for T v = 53 MeV protons. Therefore, 
proton beams are obtained with longitudinal (L) and sideways (S) polarizations at the exit 
of SOL2 by using the SOLI and SOL2 magnets, respectively. 

The beam polarization was continuously monitored by two sets of beam line polarime- 
ters (BLP1 and BLP2) Q Q after the Ring cyclotron. Each polarimeter consists of four 
conjugate-angle pairs of plastic scintillation detectors and determines the beam polarization 
via 1 H(p,p) 1 H scattering in the normal (N) and sideways (S) directions. A self-supporting 
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CH2 target with a thickness of 1.1 mg/cm 2 was used as the hydrogen target, and the elasti- 
cally scattered and recoil protons were detected in coincidence with a conjugate-angle pair 
of plastic scintillators. BLP1 and BLP2 are separated by a bending angle of 98°, and thus 
the spin precession angle in the bending magnet is about 231.1° for T p = 296 MeV protons. 
Therefore, all components (p$, Pn, Pl) of the beam polarization can be simultaneously de- 
termined using BLP1 and BLP2. The typical magnitude of the beam polarization was about 
0.70. 



B. Target and NTOF facility 

The NTOF facility [28], as illustrated in Fig. [lj consists of a beam swinger magnet, a 
neutron spin rotation (NSR) magnet, and a 100 m Time-Of-Flight (TOF) tunnel. The proton 
beam bombarded a self-supporting nat C (98.9% 12 C) target with a thickness of 89 mg/cm 2 
in the beam swinger magnet. The target energy loss was estimated to be about 270 keV by 
using the stopping power of carbon for 296 MeV protons. Neutrons from the target entered 
the TOF tunnel and were detected using NPOL3 at the end of its flight path. Protons 
downstream of the target were swept up by the beam swinger magnet and stopped by a 
graphite beam stop (Faraday cup) from which the integrated beam current was measured. 
Typical beam currents were 50 and 500 nA for the cross section and polarization transfer 
measurements, respectively. The reaction angle was changed by repositioning the target 
along the beam trajectory inside the beam swinger magnet. 

The NSR magnet was positioned at the entrance of the TOF tunnel. This magnet was 
used to precess the neutron polarization from the longitudinal V to the normal N' directions, 
enabling the measurement of the longitudinal component of the neutron polarization with 
NPOL3 as the normal component. This magnet was also used for the measurement of the 
induced polarization P [23j. In this case, the neutron polarization was precessed in the N' 
- U plane by about 120°, depending on the neutron kinetic energy. 



C. Neutron detector /polarimeter NPOL3 



The NPOL3 system 29j, illustrated in Fig. [2j consists of three planes of scintillation 
detectors. Each of the first two planes (HD1 and HD2) consists of 10 sets of one-dimensional 
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position- sensitive plastic scintillation counters (hodoscopes) with a size of 100 x 10 x 5 cm 3 , 
covering an area of 100 x 100 cm 2 . The last plane (NC) is a two-dimensional position-sensitive 
liquid scintillation counter with a size of 100 x 100 x 10 cm 3 . Both HD1 and HD2 served as 
neutron detectors and neutron polarization analyzers for the cross section and polarization 
transfer measurements, respectively, and NC acted as a catcher for the particles scattered 
by HD1 or HD2 in the polarization transfer measurements. Thin plastic scintillator planes 
(CPV and CPD) in front of HD1 and NC were used to veto and identify charged particles, 
respectively. 

The neutron energy was determined by the TOF to a given hodoscope with respect to the 
rf signal of the AVF cyclotron. A peak from l2 C(p, n) 12 N(g.s.) provided a time reference for 
the absolute timing calibration. The overall energy resolution in full width at half maximum 
(FWHM) was about 500 keV, mainly due to the target energy loss and the beam pulse width 
with contributions of about 270 and 350 keV, respectively. 

III. DATA REDUCTION 
A. Polarization observables 

A complete set of polarization observables, A y , P, and D^i = S',N',L';j = S,N,L), 
for a parity conserving reaction relates the outgoing neutron polarization p' = {jp's'iP'n'iP'l 1 ) 
to the incident proton polarization p = (ps,Pn,Pl) according to 
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The directions of the coordinate system (sideways S, normal N, and longitudinal L) are 
defined in terms of the incident proton momentum fci a b and the outgoing neutron momentum 
fc'iab in the laboratory frame as L = fc lab , V = k' lah , N = N' = (fe lab x k'i ah )/\k lah x fc' lab |, 
S = N x L, and S' = N' x L' . 
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The analyzing power A y , the induced polarization P, and the polarization transfer observ- 
able D tvjv were measured for a normally (N) polarized proton beam. The other polarization 
transfer observables, D S 's, D L i S , D s >l, and D L / L , were obtained from measurements with 
two kinds of proton beams polarized in the horizontal plane. Note that the polarization axes 
of these beams were almost orthogonal to each other. Therefore, the efficiency of measuring 
Dij is almost the same as that for pure sideways (S) and longitudinal (L) polarized proton 



beams 



23|]. 



B. Neutron detection efficiency 

The differential cross section (da/dfl) i a b in the laboratory frame is related to the observed 
neutron yield N n as 

( = N n (2) 

\dn) lah N pP AneTf h J 1 } 

where N p is the number of incident protons, p is the target density, AQ and e are the solid 
angle and intrinsic neutron detection efficiency of NPOL3 (HD1 and HD2), respectively, T 
is the neutron transmission factor along the flight path in the air, and /ii ve is the live time 
ratio. 

The product eT was determined by measuring the neutron yield from the 0° 
7 Li(jo, n) 7 Be(g.s. + 0.43 MeV) reaction which has a constant center of mass (cm.) cross sec- 



tion of (da/dfl) c , m , = 27.0±0.8 mb/sr at an incident energy range of T p = 80-795 MeV [39 ]. 
A self-supporting 7 Li (99.97%) target with a thickness of 54 mg/cm 2 was used. The eT value 
was 0.051 ± 0.003 where the uncertainty comes mainly from the uncertainties in the cross 
section for the 0° 7 Li(p, n) 7 Be(g.s. + 0.43 MeV) reaction (3%) and in the thickness of the 7 Li 
target (3%). We note that the transmission factor T has been assumed to be independent 
of reaction angle because the dependence of the flight length on reaction angle is very small 
(< 1%). 



C. Effective analyzing power 

The neutron polarization was analyzed using n + p and quasi-elastic n + C scattering 
in analyzer HD1 or HD2, and the recoiled protons were detected with catcher NC. These 
events were kinematically resolved from background events such as wrap-around and the 
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target gamma rays by using time and position information from the analyzer and catcher 
planes. Both the normal N' and sideways S' components of the neutron polarization were 
measured simultaneously, with an azimuthal distribution of n + p and quasi-elastic n + C 
events. 

The effective analyzing power A y . e ff of NPOL3 was measured using polarized neutrons 
from the GT 12 C(p, n) 12 N(g.s.) reaction at #i a b = 0°. Two kinds of polarized protons hav- 
ing normal (pn) and longitudinal (pl) polarizations were used. The corresponding neutron 
polarizations are p' N = pjyD NN (0°) and p' L = PlDll(®°) where D NN (0°) and D LL (0°) rep- 
resent the polarization transfer observables at #i a b = 0°. Then the asymmetries, and e^, 
measured by NPOL3 are 

e N = p f N^-y;e{i = PnD NN(Q°)A y ]C fi, 

\ & ) 

e L = p' L A y . cS = p L D LL (0°)A y . eS . 

As described in Sec III B\ the asymmetry was measured as the normal component using 
the NSR magnet. Because the polarization transfer observables Du(0°) for the GT transition 



satisfy [40 1 



2D NN (0°)+D LL (p°) = -l, (4) 



A,;e ff = -( 2 - + -) ( 5 ) 



A y . cS can be expressed as 

Ay.eS — — I ^ 

Pn Pl, 

using Eqs. ([3]) and (j4j). Thus, one obtains the A y . eS value without knowing ahead of time 
the Dn(0°) values. The resulting A y . e s is 0.191 ± 0.016 where the uncertainty includes the 
statistical (~ 6%) and systematic (~ 2%) uncertainties. The systematic uncertainty is 



estimated by considering the uncertainty of the beam polarization 



30|. 



D. Peak fitting 

Figure [3] shows typical excitation energy spectra of 12 C(p, n) 12 N at four momentum trans- 
fers of q = 0.14, 0.7, 1.2, and 1.7 fin -1 . The GT 1+ state at E x = MeV (ground state 
of 12 N) gives rise to a prominent peak at small momentum transfers. At large momentum 
transfers, on the other hand, its peak is small and not fully resolved from a large peak 
consisting of excited states with J n = 2 + and 2~ at E x = 0.96 and 1.19 MeV, respectively. 
Therefore, peak fitting was performed to extract the yield of the GT 1 + state. The spectra 
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were fitted at E x < 1.5 MeV where the excited states at E x = 0.96 and 1.19 MeV were 
treated as a single peak because the present energy resolution could not resolve them. The 
continuum background from wrap-around and 13 C(p, n) events was considered to be a lin- 
ear function of E x . The dashed curves in Fig. [3] show the fitting results for the individual 
peaks, while the solid curves show the sum of these contributions including the background 
indicated as the straight dotted lines. The peak fittings at all momentum transfers sufficed 
to extract the GT 1 + yield. 

IV. RESULTS AND ANALYSIS 

A. Cross section and polarization observables 

Figure H] shows the cross section for the 12 C(p, n) 12 N(g.s., 1 + ) reaction at T p = 296 MeV 
as a function of the momentum transfer q. The corresponding reaction angle # c . m . is also 
shown on the top of the figure. The momentum transfer resolution is about 0.04 fm _1 which 
is mainly due to the finite solid angle of the detector. As seen in Fig. [3X the GT state is not 
clearly resolved from the neighboring states, and thus there is a correlation between the yields 
of these two components in the peak fitting. By considering the uncertainties of the GT yields 
due to this correlation, we have estimated the systematic uncertainties. The shaded boxes 
show the uncertainties including statistical and systematic uncertainties. The statistical and 
systematic uncertainties at large momentum transfers of q ~ 1.6 fm _1 are about 2% and 
4%, respectively. A 6% uncertainty due to the cross section normalization (see Sec. IIII Bj) 
is not included. The open circles and open triangles, respectively, are data for the same 



reaction [4l| and the analogous 12 C(p,p') 12 C*(l + , T = 1) scattering [42| at T p = 295 MeV. 
The 12 C(p,p') data have been multiplied by a factor of two because of the difference in the 
isospin Clebsch-Gordan (CG) coefficients between 12 C(p, n) and 12 C(p, p'). Our data are 
consistent with the previous data within the statistical and systematic uncertainties. 

In Fig. [5], a complete set of polarization observables, D^, A y , and P, are presented for the 
12 C(p, n) 12 N(g.s., 1 + ) reaction at T p = 296 MeV as a function of momentum transfer. In the 
top right panel, the A y and P data are shown as filled and open circles, respectively, and the 
P data are offset by a momentum transfer of 0.05 fm -1 for clarity. The error bars represent 
statistical uncertainties only, while the shaded and open boxes include the systematic uncer- 
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tainties. The statistical and systematic uncertainties in for large momentum transfers 
of q > 1.0 fm -1 are about 0.19-0.24 and 0.21-0.26, respectively, which are satisfactory for 
discussing nuclear correlation effects in this momentum transfer range. 



B. DWIA calculations with shell-model wave function 



We performed microscopic DWIA calculations using the computer code DW81 [43], which 
treats the knock-on exchange amplitude exactly. Distorted waves were generated using a 



global optical mode 



T p = 20-1040 MeV 



potential (OMP) optimized for C in the proton energy range of 



44 



451 ]. with the Coulomb term turned off for the exit channel. The 



nucleon-nucleon (NN) t-matrix parameterized by Franey and Love (FL) 27| at 325 MeV 
was used as the interaction between the incident and struck nucleons. The one-body density 
matrix elements (OBDMEs) were obtained from shell-model calculations using the computer 



code oxb ash 46J . These calculations were performed in the Ohu p-shell model space using 
the Cohen-Kurath wave functions (CKWFs) [^tJ based on the (6-16) 2BME interaction. 
The transition form factor was normalized to reproduce the observed beta-decay ft value 



of 13178 s 



48] which corresponds to a GT strength B(GT) of 0.873 49j. The radial part 



of the single-particle wave functions was generated by a Woods-Saxon (WS) potential with 
r = 1.27 fm, ao = 0.67 fm [50], and a spin-orbit potential depth of V^ = 6.5 MeV j^l]. The 
depths of the WS potential were adjusted to reproduce the separation energies for the 0p 3 / 2 
orbits. 

The solid curves in Figs. 0] and [5] show the results of the calculations. The normalization 
factor for the transition form factor is N = 0.94. These calculations reproduce the experi- 
mental data reasonably well at small momentum transfers of q < 0.5 fm -1 but show poor 
agreement with the data at q > 0.5 fm" 1 . In particular, the calculations shift the momentum 
transfer dependence of the cross section to larger momentum transfers and underestimate the 
cross section at q ~ 1.6 fm - . In order to investigate the reason for this discrepancy, we next 
separated the cross section into polarized cross sections using the polarization observables. 
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C. Polarized cross sections 



The cross section / {{da/d£l) c . m . in Fig. 0J can be separated into four polarized cross 
sections IDi as 

I = ID + ID q + ID n + ID P , (6) 



where Di are the cm. polarization observables introduced by Bleszynski et al. 52j. The 
cm. coordinate system (q, n,p) is defined as q = {k' — k) /{\k' — k\), ft = (fe X k')/(\k X fc'|), 
and p = q x n, where k and k' are the momenta of the incident and outgoing nucleons 
in the cm. frame, respectively. The Di values are related to D^ in the laboratory frame 



according to 



53 1 



D = -[1 + D N N + (D S 's + D UL ) cos at 
+ {D L 's - Ds'l) sinai], 



D r 



D n 



■[1 + D NN - (D s >s + Dul) cosai 



■{D us - D S 'l) sin a 



i ) 



4 



[1 - D NN + (D S 's ~ Dl'l) cosa 2 



-{Dus + D S 'l) sina 2 ] 



(?) 



Dp = -[1 - D NN - {D S 's ~ D l >l) cosa 2 

+ {Dus + D s > L )sma 2 ], 

where a± = #i a b+fi and a 2 = 26* p — 6*i a b — fi. Here # p is the angle between k and p, and f2 is the 
relativistic spin rotation angle defined in Ref . [53] . For a plane- wave impulse approximation 
with eikonal approximation, the polarized cross sections ID { can be expressed as [53| 



ID = 


AKN D 


(|A| 2 it!o + |C| 2 itg, 


ID n = 


AKN D 


(\B\ 2 R n + \C\ 2 R ), 


ID q = 


AKN D 


{\E\ 2 R q + iD^Rp), 


ID P = 


AKN D 


(\F\ 2 R p +\D\ 2 R q ), 



(8) 
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where K is a kinematical factor, Nd is a distortion factor, A—F are the components of 
the NN t-matrix, and Ri are the response functions. Figure [6] shows the squared t-matrix 
components corresponding to each IDi. These components are derived from the FL t-matrix 
at 325 MeV. The effect of the relativistic spin rotation is so small that the D term can be 
neglected. Thus, polarized cross sections ID q and ID p represent spin-longitudinal (R g ) 
and spin-transverse (Rp) components exclusively. At forward angles, where the spin-orbit 
component \C\ 2 is very small, polarized cross sections IDq and ID n represent nonspin (Rq) 
and spin-transverse (R n ) components, respectively. 

Figure [7J shows four polarized cross sections IDi as a function of momentum transfer. 
The meaning of the error bars and shaded boxes is the same as those in Figs. H] and 
Although the present GT transition does not have a nonspin response function R , the 
nonspin polarized cross section IDq has a nonzero value due to the spin-orbit component 
\C\ 2 in Eq. (jSJ). Since ID Q is small, we will only discuss the spin-dependent polarized cross 
sections ID q , ID P , and ID n . The oscillatory pattern for the spin-longitudinal cross section 
ID q is different from those for the spin-transverse cross sections, ID P and ID n . As seen in 
Fig. EJ these patterns reflect the momentum transfer dependences of the corresponding NN t- 
matrix components. Compared with the spin-transverse t-matrix components \B\ 2 and \F\ 2 , 
the spin- longitudinal component \E\ 2 has the first minimum at lower momentum transfer of 
q ~ 0.7 fm~ . This is because the real part of E crosses zero near this momentum transfer 
due to the smallness of the pion mass. Thus the corresponding spin-longitudinal ID q shows 
the first minimum at q ~ 0.6 fm _1 . Therefore, the data verify suitable separation in the 



spin-longitudinal and spin-transverse modes based on the reaction mechanism of Ref. 53] 



The solid curves in Fig. [7J present the DWIA results with a shell-model wave function 
where the input parameters are same as those described in Sec IIVBI The calculations 
underestimate all three spin-dependent IDi at g ~ 1.6 fm _1 , and the discrepancy in the 
momentum transfer dependence is evident in the spin-transverse ID P and ID n . Thus, in 
Fig. [HI the sensitivity of the DWIA calculations for three spin-dependent IDi was investi- 
gated for changes in the parameters. Note that the solid curves are the same as those in 
Fig. [7J First, the OMP dependence of the calculations was examined by using other OMPs; 
three global OMPs (EDAD Fit 1-3) parameterized for 12 C- 208 Pb in the proton energy range 



of T p = 20-1040 MeV 



44 



45j and the OMP obtained from proton elastic scattering data on 



12 



C at T p = 318 MeV whose parameters are listed in Table HI 54| . The radial dependences of 



12 



Potential 



Vi (MeV) n (fm) a { (fm) 



Real central (i = R) 


-5.005 


1.272 


0.411 


Imaginary cetral (i = I) 


-22.55 


1.083 


0.474 


Real spin-orbit (i = RSO) 


-1.77 


0.910 


0.867 


Imaginary spin-orbit (i = ISO) 


2.71 


0.909 


0.467 


Coulomb (i = C) 




1.24 





TABLE I: The optical model parameters obtained from proton elastic scattering data on 12 C 
at T p = 318 MeV Q. The potential is defined by U(r) = V c {r) + V K f R (r) + iVifj(r) + 
[n/(m 7r c)] 2 (l/r)[VRso(d/^)/RSo(r) + iViso(d/dr)f iso (r)]((T x L), where V c is the coulomb po- 
tential for a uniformly charged sphere and f%(r) = [1 + exp{(r — r^ 1 / 3 )/^}] -1 . 

these OMPs for the incident channel are shown in Fig. [HJ and the DWIA results are shown 
in Fig. [8](a) by the bands. The OMP dependence of the spin-longitudinal ID q is small, 
whereas those of the spin-transverse ID p and ID n are significantly larger near the cross 
section minimum at q ~ 1.4 fm -1 . We also performed DWIA calculations using the neutron 
global OMPs for 12 C- 238 U in the neutron energy range of T n = 20-1000 MeV (s5 1 for the 
exit channel, and the results are plotted as the dashed curves in Fig. E](a). The use of the 
neutron global OMPs gives larger values near the cross section minimum at q > 1.4 fm -1 in 
the spin-transverse mode. However, neither the discrepancy in the angular distribution nor 
the underestimation in the cross section at large momentum transfers can be explained by 
the OMP uncertainties. 

Second, the proton-particle and neutron-hole configuration dependences were investi- 
gated. The bands in Fig. MJo) are DWIA results with other CKWFs based on the (8-16) 
2BME and (8-16) POT interactions 0|. We also performed DWIA calculations for a pure 
OP1/2OP3/2 transition from the Hartree-Fock (HF) state of 12 C (the state fully occupying the 
Osi/2 and 0p 3 / 2 orbits), and the results are shown as the dashed curves. Table ILT1 summarizes 
the OBDMEs and -B(GT) together with the corresponding normalization factors for the 
transition form factors. The configuration dependence is small for all three IDi, and thus the 
discrepancy between the experimental and theoretical results is not resolved by considering 
the configuration dependence. 
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shell-model wave function 




OBDME 




5(GT) 


N 




0pi/ 2 0pi L 


0»i /o0p„ } n 


OP^/oOPt Ir, 


0^3/20^0/!, 






(6~16)2BME 


0.0859 


0.6672 


0.3228 


0.0925 


0.929 


0.94 


(8-16)2BME 


0.0733 


0.6915 


0.3262 


0.0822 


0.992 


0.88 


(8-16)POT 


0.0582 


0.6902 


0.3393 


0.0763 


0.921 


0.95 


Pure Opx^OpsJz 


0.0 


1.0 


0.0 


0.0 


5.228 


0.17 



TABLE II: One-body density matrix elements and Gamow- Teller strengths for the 
12 C(p, n) 12 N(g.s., 1 + ) reaction used in the DWIA calculations. The normalization factors for the 
transition form factors are also listed. 

Finally, we investigated the dependence on the radial wave function. The dashed curves in 
Fig. [S^c) are DWIA results using a harmonic oscillator (HO) potential with a size parameter 
of b = 1.53 fm _1 . This parameter was obtained from an analysis of the electron scattering 



on 12 C to the stretched 4~, T = 1 state at E x = 19.55 MeV 56( with the center-of-mass 
correction taken into account [57J]. The results are almost the same as those for the WS 
potential. Thus the discrepancy between the experimental and theoretical results could not 
be explained by the radial wave function dependence. 

Figure (TU] represents the parameter dependence for the orthogonal components of the 
polarization transfer observables, Dnn (left panels), Ds>s (middle panels), and Dl>l (right 
panels). Based on the calculations in Figs. [8] and [101 the experimental data at large mo- 
mentum transfers cannot be reproduced within the framework of the DWIA employing 
shell-model wave functions. Therefore, in the following section, nuclear correlation effects 
beyond the shell model are investigated. 

V. DISCUSSION 

In this section, the experimental spin-dependent polarized cross sections are compared 
with the DWIA calculations using RPA response functions in order to investigate the nuclear 
correlation effects beyond the shell model. 
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A. DWIA+RPA calculations 



We performed DWIA+RPA calculations using the computer code CRDW 58] . The for- 



malism of the calculations is discussed in Refs. [51|, [58j. The spin-isospin response functions 
were calculated in a continuum RPA including the A degrees of freedom. We further uti- 
lized a ring approximation Jssj], and used the tc + p + g' model interaction for the effective 
interaction, which is expressed as [l] 



V e g(q, u) = V L (q, u) + V T {q, u), 



(9) 



where Vl and Vr are the spin-longitudinal and spin-transverse effective interactions, re- 
spectively. They are determined by the pion and rho-meson exchange interactions and the 
Landau-Migdal (LM) interaction specified by the LM parameters, g' NN , g'^^, and g'^, as 



V L (q,u) 



f 2 f a 2 



2 \ 9'nn + , l2 _*2_ rn 2 T Inn(Q^) ) ( CT i • 0)^2 ■ q)(r 1 ■ T 2 ) 



— — 



. fnNNfirNA ( / Q n / s-p / 

H 5 9na H 5 o T 1 ttNnKQ, ^v 1 nNA\Q) ^ , 

m z y _ gZ _ m z 

x [{(<n ■ q)(S 2 ■ q)(ri • T 2 ) + (1 <-> 2)} + h.c] 



f2 / „2 

J 7T 



m 



Q 



X 



• I S'aa "I 2 2 ^ijva(?) w 
~ \ or — q z — m z 

5x • $)(Sj • q)(T! • T\) + (5x • <z)(S 2 • q)(T! • T 2 ; 



and 

V T (q,uj) 



ft 



NN I i 



mi 



9nn + 



x(cr 1 x q)(<x 2 x q)(n • t 2 ) 

, firNNfwNA ( , n 
^ n I <?iVA "+" W - ^" 



lu^ — q 2 — m 2 



x [{(<r i x q)(5 2 x q){ Tl ■ T 2 ) + (1 <-> 2)} + h.c] 



~2 ( ^AA + Cp, ,2 _ „2 _ ^2 r pA f A(?' W , 



ni: 



h.c. 



or — g — m p 

5 X x g)(flft x q)( Tl . T t) + ( 5l x ^(^ x g)(T x • T 2 )| + h.c. 



(10) 



(11) 
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where and m p are the pion and rho-meson masses, a and r are the spin and isospin 
operators of the nucleon N, and S and T are the spin and isospin transition operators 
from iV to A. The coupling constants and meson parameters for the pion and rho-meson 



exchange interactions from a Bonn potential were used, which treats A explicitly [59|. The 
LM interaction effectively represents the short range correlations and the exchange terms 
in the RPA, and t he L M parameters have been estimated to be g' N N = 0.65 ± 0.15 and 
9na = 0-35 ± 0.15 [60] by using thepeak position of the GT giant resonance and the GT 
quenching factor at q = fm _1 2], [J, as well as the isovector spin-longitudinal polarized 
cross section in the QES process at q ~ 1.7 fm -1 [25]. We fixed g' AA = 0.5 61[ because the 
g' AA dependence in the results is weak. The response functions are normalized to reproduce 
the experimental -B(GT). 

The ground state of 12 C was assumed to be a HF state. However, as seen in Table HT| the 
ground state correlation which is included in the shell-model calculations plays an important 
role in reproducing the experimental B(GT) value of 0.873. In order to include the shell- 
model (configuration-mixing) effects effectively, we used much smaller normalization factors 
iV than those of the shell-model calculations, namely N = 0.28 and N = 0.17 in the 
calculations with and without RPA correlations, respectively. 

The nonlocality of the nuclear mean field was treated using the local effective mass 
approximation [1] in the form 

m*(r) =m s - ^U(™n ~ m*(0)), (12) 
where is the nucleon mass and fws{ r ) is a WS radial form. Here we adapted the standard 



value of m*(0) = 0.7rriN [62j, [63] 



We used the same OMPs and single-particle wave functions as described in Sec. IIVBI 
The iViV t-matrix parameterized by Franey and Love at 325 MeV was used, and the ex- 
change terms were approximated by contact terms following the prescription by Love and 



Franey 



64]. 



In Fig. HU the consistency between theoretical results using the computer codes CRDW 
and DW81 is checked. The solid curves are the CRDW results with a free response function 
employing m*(0) = tun, whereas the dashed curves are the DW81 results with the corre- 
sponding wave function as for the dashed curves in Fig. MJo). The normalization factor iV 
of 0.17 is common, and both calculations for the unpolarized cross section / are in good 
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agreement with each other. The theoretical results for IDi show good consistency except 
near the cross section minimum at q ~ 1.4 fm _1 . Therefore, in the next subsection, the 
experimental data are compared with the DWIA calculations using CRDW in order to inves- 
tigate nuclear structure effects which could not be included in the preceding calculations, 
such as RPA correlations. 



B. Comparison with DWIA+RPA calculations 

First, the nonlocality of the nuclear mean field was investigated. The dotted and dashed 
curves in Fig. [T2l show the DWIA results with a free response function using m*(0) = rriN 
and m*(0) = 0.7rriN, respectively, and N = 0.17. The angular distributions of all three spin- 
dependent IDi curves shift to lower momentum transfer due to the nonlocality of the nuclear 
mean field, so that the agreement with the data is improved. This shift arises because the 



transition form factor moves outward due to the Perey effect 65(. However, there remains 



a discrepancy between the experimental and theoretical results at around q ~ 1.6 fm _1 . 

Next we considered the nuclear correlation effects in the RPA. The solid curves in Fig. 
show the results of DWIA+RPA calculations using g' NN = 0.65, g' NA = 0.35, and m*(0) = 
0.7mjv with N = 0.28. The bands represent the g' NN and g' NA dependences with g' NN = 
0.65±0.15 and g' NA = 0.35±0.15. In the continuum RPA, the GT state couples to particle- 
unbound 1 + states, which shifts the response function in coordinate space to larger r values. 
Thus the angular distributions further shift to lower momentum transfer. Furthermore, 
the RPA correlation enhances all three IDi at large momentum transfers of q ~ 1.6 fm~ , 
improving the agreement with the data. In the analysis of QES, the spin-transverse ID P 



and ID n are quenched due to the repulsion of the spin-transverse interaction Vr [231 . 125 ] . 
However, in the analysis shown in Fig. [T21 the transition form factors are normalized to 
reproduce -B(GT). This means that the quenching due to the repulsive LM interaction is 



effectively included through the normalization factor N [32J. Therefore, in Fig. [12], the 
attractive rho-meson exchange effects are seen as an enhancement of ID p and ID n at q ~ 
1.6 fm _1 Also note that the modification of the momentum transfer dependences (due to 
the shape change of the response functions in r-space which is not included through the 
normalization) is important for the magnitude of IDi. 



In some theoretical studies 10, IfiJa, [67j, the LM parameters are taken to have the 
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momentum- and energy-transfer dependence. We investigated this effect on IDi by us- 
ing the effective interaction V e s by Alberico et al. [10(j , in which the LM parameters have the 
momentum- and energy-transfer dependence with the dipole form factors. The DWIA+RPA 
results with this effective interaction are shown in Fig. [T3J The dashed and solid curves cor- 
respond to the calculations with and without the dipole form factors, respectively. We 
note that the coupling constants and meson parameters are slightly different from those 



in Ref. 59j, and thus the results without the dipole form factors are also slightly different 
from those shown in Fig. [121 For the spin-longitudinal mode, the use of the dipole form 
factors enhances ID q at large momentum transfers of q ~ 1.6 fm _1 , which is due to the 
more attractive spin-longitudinal interaction Vj,. However, the effect is small, and the same 
results can be achieved by using smaller and reasonable LM parameters of g' NN ~ 0.55 and 
g' NA ~ 0.30. For the spin-transverse mode, the form factor effects are very small in both 
I D p and ID n because of the insensitivity to the spin-transverse interaction Vr |l|. 

Figure [14] compares the experimental and theoretical results for the cross section and or- 
thogonal components of the polarization transfer observables, Dnn, Ds>s, and Dl>l- These 
quantities are better reproduced by considering RPA correlations together with the nonlo- 
cality of the nuclear mean field, particularly for the momentum transfer dependences. From 
the analyses in Figs. [12] and [HI we conclude that our data support the existence of pionic 
and rho-mesonic correlations in nuclei at large momentum transfers. 

VI. FINAL REMARKS 

The DWIA calculations including RPA correlations reproduce the experimental data for 
the spin-longitudinal cross section ID q and give improved descriptions of the spin-transverse 
cross sections, ID P and ID n . However, the experimental values of ID P and ID n remain 
larger than the calculated values at momentum transfers of q ~ 1.6 fm" 1 by factors of about 
1.4 and 2.0, respectively. The magnitudes of the observed enhancements are significantly 
different for ID p and ID n . Since the spin-transverse response is common for ID p and ID n , 
medium modifications of the effective NN interaction are considered as a possible reason 
for the observed enhancements. Such modifications have been discussed using data for the 



stretched state excitations in Refs. 



31 



68, 



69], where it is reported that the experimental 



values of ID n at large momentum transfers are increased in magnitude by a factor of about 
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1.5, indicating an enhancement in the corresponding NN scattering amplitude B. Thus, a 
larger value of B might be responsible for the enhancement in the experimental values of 
ID n from Eq. (181) . However, since no modification of the NN amplitude F is observed in the 
stretched state excitations, the enhancement in ID p cannot be explained by modifications 
of the NN interaction. We note that medium modifications of the effective NN interaction 
have been also discussed in the relativistic framework [3]. If one takes into account the 
increase in B, the enhancement of ID n is reduced to 1.3, which is almost same as that of 
ID p . Consequently the experimental values of ID p and ID n may be en 
mechanism, such as higher order (e.g., 2p-2h) configuration mixing 
comprehensive and detailed theoretical analyses are needed including effects of higher order 
configuration mixing and medium modifications of the effective NN interaction. 



lanced by a common 

n n 

71, 72| . Thus more 



VII. SUMMARY AND CONCLUSION 



We have measured differential cross sections and complete sets of polarization observables 
for the 12 C(p, n) 12 N(g.s.,l + ) reaction at T p = 296 MeV with momentum transfers of q = 0.1 
to 2.2 fm _1 in order to investigate nuclear correlation effects inside the nuclei. The experi- 
mental polarized cross sections IDi have been compared with DWIA calculations employing 
shell-model wave functions. For all three spin-dependent IDi, a significant difference in the 
momentum-transfer dependence and an enhancement around q ~ 1.6 fm _1 were observed 
when compared to calculations. The use of a local effective mass of m*(0) = 0.7m^ im- 
proves the agreement with the data, but still underestimates the cross section at around 
q ~ 1.6 fm . These underestimations for all three spin-dependent IDi are partly resolved 
by DWIA calculations employing an RPA response function with g' NN = 0.65, g' NA = 0.35, 
and m*(0) = 0.7m7v, supporting the existence of pionic and rho-mesonic correlations in the 
nuclei. This finding is the first indication for observing pionic and rho-mesonic correlation 
effects separately. To understand the nuclear correlation effects quantitatively, theoretical 
analyses are required that include effects of higher order configuration mixing and medium 
modifications in the effective NN interaction. 
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Beam 



FIG. 1: A schematic layout of the NTOF facility (not to scale). The coordinate systems for incident 
protons and outgoing neutrons are also shown. S (Sideways), N (Normal) and L (Longitudinal) 
form a right-handed system for incident protons and S', N' , L' for outgoing neutrons. 
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FIG. 2: (Color online) A schematic view of the neutron detector /polarimeter NPOL3. In the 
detector mode, HD1 and HD2 act as neutron detectors. In the polarimetry mode, HD1 and HD2 
are the analyzer planes while NC is the catcher plane. Thin plastic scintillator planes are used to 
veto (CPV) or identify (CPD) charged particles. 
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FIG. 3: (Color online) Excitation energy spectra for 12 C(p, n) 12 N at T p = 296 MeV and q = 0.14, 
0.7, 1.2, and 1.7 far 1 . The dashed curves are fits to the individual peaks. The solid curves indicate 
the sum of the peak contributions including the background plotted as the straight dotted lines. 
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FIG. 4: (Color online) Measured cross sections for 12 C(p, n) 12 N(g.s., 1 + ) at T p = 296 MeV (filled 
circles) as a function of momentum transfer. The corresponding reaction angle # c . m . is also shown 
on the top of the figure. The open circles are data at T p = 295 MeV [41]. The open triangles are 



12, 



C(p,p') 12 C*(1+,T= 1) data at T p = 295 MeV 



421 ]. multiplied by a factor of two as described in 



the text. The solid curve shows DWIA calculations using a shell-model wave function. 
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FIG. 5: (Color online) Measured polarization transfer observables Dij, analyzing power A y , and 
induced polarization P for 12 C(p, n) 12 N(g.s., 1 + ) at T p = 296 MeV. The induced polarization data 
P are offset by a momentum transfer of 0.05 fm _1 so that the analyzing power A y and induced 
polarization P data can be distinguished. The solid and dashed curves are the results of DWIA 
calculations with a shell-model wave function. 
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FIG. 6: (Color online) Squared £-matrix components calculated from the FL i-matrix at 325 MeV. 
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FIG. 7: (Color online) Measured polarized cross sections for 12 C(p, n) 12 N(g.s., 1+) at T p = 
296 MeV. The solid curves are the results of DWIA calculations using the shell-model wave 
function. 
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FIG. 8: (Color online) Parameter dependence of the calculations for the spin-dependent polarized 
cross sections ID q (left panels), ID P (middle panels), and ID n (right panels). The solid curves 
are the same as those in Fig. [71 The bands and dashed curves present DWIA results with other 
parameters: (a) four different OMP parameters (bands) and neutron global OMPs for the exit 
channel (dashed curves); (b) two different CKWFs (bands) and a pure Opx^Op^ configuration 
(dashed curves); and (c) the HO potential (dashed curves). 
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FIG. 9: (Color online) The radial dependences of the OMPs for the incident channel used in 
the DWIA calculations. The solid curves and the bands show the global OMPs optimized for 
12 C (EDM) and 12 C- 208 Pb (EDM Fit 1-3) in the proton energy range of T p = 20-1040 MeV, 



respectively 



44 



451 ]. The dashed curves represent the OMP obtained from proton elastic scattering 



data on 12 C at T p = 318 MeV 



54J. 



33 



1.0 
0.5 
0.0 
-0.5 
-1.0 

1.0 
0.5 
0.0 
-0.5 
-1.0 

1.0 
0.5 
0.0 
-0.5 
-1.0 



■ (a) D> 

Jft, ( 


IN \ 
> 




■Dl>l 

1 




.... 1 .... 1 .... I .... T ... 


.... 1 .... 1 .... 1 .... I ... 


-v 

.... 1 .... 1 .. . 


A < 

f\\l 
1 \ 




:(b) 


► 


,,,,,,,,, 1 ,,,,,,,,,,,, , 

< '' « A I : 


,,,,,,,,, 1 ,, , 

1 

: -t-« 




■ *m 


: \j/ : 


^^^^^^ 






"(c) 


' ' ' . 

> 


■ , , 1 1 1 1 1 , 1 1 1 1 1 1 , , i , 1 1 1 , 


■ , , 1 1 1 1 1 , 1 1 1 1 1 1 , , i , 1 1 1 , 






.... 1 .... 1 .. . 





0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

Momentum transfer q (fm -1 ) 



FIG. 10: (Color online) Same as in Fig. [8] but for the diagonal polarization transfer observables: 
Dnn (left panels), D$>s (middle panels), and Dl>l (right panels). 
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FIG. 11: (Color online) Comparison between calculations using the computer codes crdw and 
dw81 for the cross section / and polarized cross sections, ID q , ID P , and ID n . The solid and dashed 
curves represent DWIA calculations using CRDW and dw81, respectively. The normalization factor 
iV of 0.17 is common for both calculations. 
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FIG. 12: (Color online) Comparison between experimental and theoretical results of polarized 
cross sections ID q , ID P , and ID n for 12 C(p, n) 12 N(g.s., 1+) at T p = 296 MeV. The dotted and 
dashed curves present the DWIA results with a free response function using m*(0) = tun and 
m*(0) = 0.7mjy, respectively, and ./V = 0.17. The solid curves denote DWIA results for an RPA 
response function with {g' NN , fl^A' 9aa) = (0.65,0.35,0.50) , m*(0) = 0.7mjv and N = 0.28. 
The bands are the g' NN and g' NA dependences of the DWIA results with g' NN = 0.65 ± 0.15 and 
g' NA = 0.35 ±0.15. 
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FIG. 13: (Color online) DWIA predictions for 12 C(p, n) 12 N(g.s., 1+) at T p = 296 MeV. The dashed 
and solid curves are the DWIA results with the RPA response functions using the LM parameters 
with and without the dipole form factors, respectively. 
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FIG. 14: (Color online) Same as in Fig. [12] but for the cross section and diagonal polarization 
transfer observables, -D/vat, As's, an d 'Dul- 
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